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ABSTRACT 

It is found that IE 1207.4-5209 could be a low-mass bare strange star if 
its small radius or low altitude cyclotron formation can be identified. The age 
problems of five sources could be solved by a fossil-disk-assisted torque. The 
magnetic dipole radiation dominates the evolution of PSR B1757-24 at present, 
and the others are in propeller (or tracking) phases. 

Subject headings: pulsars: general — pulsars: individual (IE 1207.4-5209) — 
stars: neutron 



1. Introduction 

The age of neutron star is an essential parameter, which is relevant to the physics of 
supernova explosion and thereafter the evolution of stars. However, it is still a big problem 
now to determine generally an exact value of age (except the crab pulsar). It is a conven- 
tional and convenient way to obtain the age for rotation-powered neutron stars by equalizing 
the energy lose rate of spindown to that of magnetodipole radiation, assuming that the in- 
clination angle between magnetic and rotational axes is a = 90° (e.g., Manchester & Taylor 
1977). The conclusion keeps quantitatively for any a, as long as the braking torques due to 
magnetodipole radiation and the unipolar generator are combined (Xu & Qiao 2001). The 
resultant age, the so-called characteristic age, is Tc = P/ {2P) if the initial period Pq is much 
smaller than the present period P. The age, T^, is generally considered as the true one for 
a neutron star with P > 100 ms since most newborn neutron stars could rotate initially at 
Po ~ (20 - 30) ms (e.g., Xu et al. 2002). 

It challenges the opinion above that the ages of a few supernova remnants are incon- 
sistent with Tc of their related isolated stars (Table 1), which implies that some additional 
torque mechanisms do contribute to star braking. Among the five stars, three of them have 
Tc > lOTsNR (the age of supernova remnant), and other two Tc < 2Tsnr- In addition, elec- 
tron cyclotron resonant lines are detected in two of the five neutron stars (IE 1207.4-5209: 
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Bignami et al. 2003; IE 2259+568: Iwasawa et al. 1992), but the inferred magnetic fields, 
-Bcyc, from which are significantly smaller than that in the magnetic dipole radiation model, 
B^. The most prominent one in this age discrepancy issue is IE 1207.4-5209, which has 
To > SOTsNR and B^y, < BJ30. 

One probable and popularly-discussed way to solve the age problem is through an ad- 
ditional accretion torque (Marsden et al. 2001, Alpar et al. 2001, Menou et al. 2001). In 
this paper, whether an additional accretion torque can possibly solve the age discrepancy is 
investigated, including discussions about possible astrophysical implications. 

2. The case of IE 1207.4-5209 

The key point in the age discrepancy of IE 1207.4-5209 is how to spin down from 
Pq ~ 20 ms to 424 ms in a short time of Tsnr ~ 7 kyrs if its true age is Tsnr. Certainly the 
problem disappears if one assumes a long initial period Pq ~ 400 ms or a large braking index 
n ~ 50 (Pavlov ct al. 2002); but this is not of Occam's razor since it is generally believed 
that rotation-powered radio pulsars born with ~ 20 ms, and brake with index < 3. May an 
additional accretion torque help the spindown? Actually, in an effort to reconcile B^ with 
Bcyc, an accretion model for IE 1207.4-5209 was proposed (Xu et al. 2003). However, a 
very difficulty in the model is how to choose a time-dependent accretion rate Md(i), and to 
determine the propeller torque with the rate M^- 

Nevertheless, the propeller phase works in the centrifugal inhibition regime when rm > 
Tc, for a star with mass M and magnetic moment the corotation radius Tc = [GM/ (47r^)]^/^P^/^ 
and the magnetospheric radius — {p? / {M^\/2GM)Y'/'^ . To avoid the complex calculations 
of magnetohydrodynamics, the rotation energy loss due to propeller torque could be simply 

introduced as E^, = —GM^M/R^, based on the energy conservation law. This is unphysical, 
but should be an limit for accretion braking. As rm decreases {-^ r^), increases, and \Eg\ 
increases too. Therefore the most efficient spindown (MESD) takes place when rm ^ r^. 

In a model where the propeller and electromagnetic torques are combined, in the MESD 
case, one can derive the period evolution 

P < l.lBl^Rt{M/MQ)-\t/jTs) + Po (ms), (1) 

where B12 = i?/(10^^G) and = it!/(10^cm). The right hand of Eq.(l) is an upper limit 
of P because 1, a realistic accretion rate may not be as high as that of MESD; and 2, the 
corresponding braking torque is not so effective. If IE 1207.4-5209 is a conventional neutron 
star with mass ~ Mq and radius ~ 10^ cm, and the line features are related to cyclotron 
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absorptions near the surface (Xu et al. 2003; the polar magnetic field is thus 6 x 10^° G), 
one has P < 3.8(i/kyrs) + Pq (ms). 

Therefore, assuming IE 1207.4-5209 has a true age t ~ 7 kyrs and an initial period 
Pq ~ 20 ms, the upper hmit of the present period is ~ 40 ms (<^ P — 434 ms), and the age 
discrepancy can then not be solved in the conventional neutron star model. However, if IE 
1207.4-5209 is a strange star with low mass, for instance R — 1 km (and the mass is thus 
~ 1O~^M0, since low-mass strange stars have almost a homogenous density ~ 4 x lO^'^g/cm^; 
Alcock et al. 1986), the upper limit is then P ^ 110Bl2{t/jTs) + Pq (ms). In this case, IE 
1207.4-5209 could spin down to ~ 2.8 s during ~ 7 kyrs if its polar magnetic field 6 x 10"^° 
G. In fact, the fitted radius of IE 1207.4-5209 with a blackbody model is only ~ 1km 
(Mereghetti et al. 1996; Vasisht et al. 1997) although a lager radius is possible if a light- 
element atmosphere is applied (Zavhn, Pavlov & Triimper 1998). The best- fit tow blackbody 
model of XMM-Newton data indicates an emitting radius ~ 3 km for the soft component 
with temperature ~ 200 eV (Bignami et al. 2003). Combined with its non-atomic feature 
spectrum, we may suggest that IE 1207.4-5209 is a low-mass strange star with bare quark 
surface (Xu 2002, Xu et al. 2003). 

An alternative possibihty is that IE 1207.4-5209 is a conventional neutron star, but 
the cyclotron resonant absorption forms far away from the surface. The polar magnetic 
field is ~ 6 X 10^°[(i? + h)/R]^G if the resonant lines form at a height h. From Eq.(l), 
424 < 1.1 X 7 X 10^-8^2 + 20, we estimate a low limit of the polar magnetic field to be 
~ 2.3 X 10^^ G. This implies that the resonant absorption region should be at a level of > 16 
km height from the surface. Certainly, in case of no propeller torque (i.e., the polar magnetic 
field is (1.7 — 3.6) x lO^^G), the height of resonant absorption region is (30 — 40)km. 

It is worth noting that IE 1207.4-5209 could be a low-mass neutron star with polar 
magnetic field B12 = 0.06. From Eq.(l) and the conditions of MESD, one has i?g > 4:{M/Mq) 
for Pq ~ 20 ms. This implies a neutron star with radius i? > 10 km but mass M < Mq (e.g., 
Shapiro & Teukolsky 1983). This result may have difficulties in explaining 1, a non-atomic 
spectrum (Xu et al. 2003), and 2, a possible small radius observed (Mereghetti et al. 1996; 
Vasisht et al. 1997; Bignami et al. 2003) and even the fitting result of a neutron star with 
10 km and 1.4M0 (Zavlin et al. 1998). 



3. Other sources 

If other sources fist in Table 1 are neutron stars with Rq — 1 and M — Mq, the low 
limits of polar magnetic fields are 6.1 x 10^^ G for IE 2259+586, 4.9 x 10^° G for PSR 
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B1757-24, 1.6 x 10^^ G for PSR J1811-1925, and (2.5 - 5.6) x 10^^ G for PSR J1846-0258. 

Among these sources, only possible cyclotron absorption is found in IE 2259+586, and the 
limit field is within the range of that inferred from cyclotron line. This suggests that the 
cyclotron resonant may take place just above the stellar surface. 

An interesting question is: how much mass could accrete during the propeller phase in 
case of MESD? Certainly, only a very small part of this matter can accrete onto the stellar 
surface. When MESD works, one obtains the accretion rate 

from Tm ~ Tc- If the quantities are re-scaled, m = M^/Mq, t = t/ji, and p = P/ms, one 
has dm/dr ~ 0.24 /x|q(M/M0)~^/^p~^/^. Combining with Eq.(l), one comes to 

poo 

m< dr = 0.16(M/MQ)-2/3j^5Po'^', (2) 

where po — -Po/ms. Note that the upper limit of accretion mass, Md, in the right hand of 
Eq.(2) does not depend on the magnetic momentum. Typically for po = 20, the upper limit 

of accretion mass is 2.95 x 1O~^M0, which is reasonable since the amount of the fall-back 
material after supernova explosion could be as high as O.IM© (Lin, Wooslcy & Bodenheimer 
1991; Chevalier 1989). Due to r-mode instability, a nascent neutron star may loss rapidly its 
angular momentum though gravitational radiation if the initial period is less than ~ 3 — 5 
ms (Andersson & Kokkotas 2001). The upper limit of accretion mass in case of MESD is 
~ O.O4M0 for po = 3. These results indicate that the fall-back matter is enough to brake 
the center stars by propeller torque in MESD case. 



4. A model with self-similar accretion rate 

Though the study about MESD torque provides some useful information on the accretion 
model, including the appropriate magnetic field and the mass of fall-back disk around a 
neutron star, the accretion rate of MESD torque is questionable in reahstic cases. After a 
dynamical time a fossil disk may form. For a viscosity-driven disk, the accretion could be in 
a self-similar way, with an accretion rate of (Cannizzo, Lee & Goodman 1990) 

m = mo, < t < T; m = mQ{t/T)-", t > T, (3) 

where T is of order the dynamical time, m = dm/dt. Assuming T ~ 1ms, an initial disk 
mass ~ O.OO6M0, and an opacity dominated by electron scattering {a — 7/6), Chatterjee, 
Hernquist and Narayan (2000) develop a first detailed model of fossil disk accretion for 
anomalous X-ray pulsars (AXPs). However it is noted by Francischelli and Wijers (2002) that 
Kramers opacity may prevail in the fossil disk (i.e., a = 1.25). In the regime of conventional 
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neutron stars, we will calculate the accretion torque through the realistic accretion rate of 
Eq.(3), assuming a = 1.25 and T = 1ms, with the inclusion of magnetic dipole radiation. 
The spinup/spindown torque proposed by Menou et al. (1999), 

j = 2My^n^{r^)[i - n/n^{r^)], (4) 

is applied for the action of fossil disk in the model, where fl — 27:/ P, Ok(?"m) is the Keplerian 
angular velocity at the magnetospheric boundary. 

One may compute the accretion rate, rfi, as well as the spin evolution P{t). The total 
disk mass could be m = rhdt. It is worth noting that the disk mass obtained in this way 
could be much larger than that in MESD case because of the inclusion of the high accretion 
in the initial period. 

We think that the accretion rate characterized by Eq.(3) is on average in a sense. The 
period derivative, P, may be affected by dynamical instabilities or some stochastic processes, 
whereas the period, P, is of the integration over a very long time. We therefore calculate 
-P(^snr) for any disk mass, m, and polar magnetic field, B12, of neutron stars. For IE 
1207.4-5209, the calculated contour of period relative error, |P(Tsnr) - P\/P {P = 0.424s 
for IE 1207.4-5209), is shown in Fig.l. A reasonable parameter set (disk mass m and polar 
magnetic field B12) is chosen if the following criteria are met: 1, the period relative error 
is smaller; 2, m < 0.1; 3, Bi2,cyc < -B12 < -Bi2,d (Table 1). We then have m — 0.054 
and B12 — 3.55. The parameter sets for other sources can also be obtained in this way, 
which are listed on Fig.2 (except for PSR J1811-1925). The period evolution curves, with 
these parameters, are drawn in Fig.2. Note these curves do not change significantly if the 
parameter sets shift reasonably. 

The heights, h, of cyclotron resonant scattering regions can be obtained, based on the 
differences of parametric magnetic field, i?i2, and the field inferred from absorption features, 
5i2,cyc. It is found that /i ~ 29 km and ~ (8.8-15) km for IE 1207.4-5209 and IE 2259+586, 
respectively, in the model. 

Whether or not the disk will influence the spindown of the neutron star or suppress the 
radio emission will depend on the location of rm relative to the light cylinder radius, tl, and 
the corotation radius (Chattcrjcc ct al. 2000). Magnetic dipole radiation dominates, and 
the disk and the star will effectively evolve independently if > tl; but in other cases, 
accretion onto the star will lead to accretion- induced X-ray emission with radio quiet. We 
see from Fig.2 that the condition of Tm > tl is satisfied only for PSR B1757-24 when it is 
older than ~ 10^ years. We are therefore not surprise that PSR B1757-24 is now radio loud 
whereas the others (IE 1207.4-5209, IE 2259+586, and PSR J1846-0258) are radio quiet. 
The AXP IE 2259+586 is in a tracking phase, and we expect other two (IE 1207.4-5209 and 
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PSR J1846-0258) will evolve to be AXPs when they are in tracking phases too. 

PSR J1811-1925 is an interesting exception among the five sources, whose age is certain 
if it has physical association with the remnant of a supernova recorded in A.D. 386. In its 
calculated contour, we can only choose {B12 — 2.6, m = 0.1} or {B12 — 1.6, m = 2.2}; both 
the parameter sets are not reasonable (i.e., can not meet those 3 criteria). This may imply 
that the accretion of PSR J1811-1925 is not self-similar. Recalling that the low limit of polar 
field is only 1.6 x 10"G if in MESD case, we could suggest that PSR J1811-1925 has a field 
within (1.6 — 17) x lO^^G, with an accretion stronger than that of Eq.(3) but weaker than 
that in MESD case. This result hints that PSR J181 1-1925 is radio quiet (Crawford et al. 
1998). In addition, the parametric field, Bu = 3.55, chosen for IE 1207.4-5209, which is 
close to -Bi2,d = 3.6, would also indicate that the real accretion is not described by Eq.(3). 
In fact the accretion of Eq.(3) is for the capture of material by black holes where magnetic 
field is not important, which could differ from that for neutron stars with strong fields. 



5. Conclusions and Discussions 

The possibility of solving the age discrepancy by an accretion-assisted torque is dis- 
cussed. We find that: 1, IE 1207.4-5209 can not be a neutron star, but a low-mass bare 
strange star, if the cyclotron resonant region is near the polar cap with a magnetic field of 
6 X lO^^G; whereas it could be a conventional neutron star if the cyclotron lines form at a 
height of (16 — 40)km. An identification of smaller radius or low altitude cyclotron formation 
favors a low-mass bare strange star model for IE 1207.4-5209. 2, Among the five sources 
with age problems, the magnetic dipole radiation dominates the evolution of PSR B1757-24 
at present, and the others are in propeller (or tracking) phases. 3, The real accretion around 
these sources may differ from a self-similar one (Eq.(3)), at least for PSR J1811-1925. 4, 
By a calculation with self-similar accretion, it is suggested that PSR J1846-0258 and IE 
1207.4-5209 (and probably PSR J1811-1925) would evolve to be anomalous X-ray pulsars in 
the future. 

The debris disks formed following supernovae are currently conjectured also for inter- 
preting other astrophysical phenomena, e.g, anomalous X-ray pulsars and soft 7— ray re- 
peaters. Factually, these disks around the sources could be bright in a wide spectral range. 
Recent discoveries of possible optical and near-infrared emission from a few AXPs may be 
a hint of such kind fallback accretion disks (IE 2259+586: Hulleman et al. 2001; IRXS 
J170849-400910: Israel et al. 2002; IE 1048.1-5937: Wang & Chakrabarty 2002). Although 
a comparison of optical and near-infrared observations with theoretical predictions of spec- 
tra of disks around neutron stars (Perna et al. 2000) have helped rule out the presence of 
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disks in some cases, more detailed studies in this aspect is still necessary, which may be an 
effective way to test the fossil-disk model for young neutron stars. 
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Table 1. List of neutron stars with age discrepancies 



stars SNRs P(s) P(10-i3s/s) Tc(103y) TsNR(103y) BJj.d 'Sf2,cyc ^^ef. 



IE 1207.4-5209 


PKS 1209-51/52 


0.424 


0.07-0.3 


200-900 


~ 7 


1.7-3.6 


0.06 


1,2 


IE 2259+586t 


CTB 109 


6.98 


4.84 


228 


17 


59 


0.4-0.9 


3,4,5 


PSR B1757-24 


G5. 4-1.2 


0.125 


1.28 


16 


> 39 


4.0 




6 


PSR .11811-1925 


Gll.2-0..'? 


0.065 


0.41 


21 


1.6 


1.7 




7 


PSR .llSl(i-02.-)S 


Kcs 7.") 


l).:-!25 


71 


l).72 


0.9- 1.:-! 


19 




s 



*Thc magnetic fields in the magnetic dipole radiation model; B,i = .'?.2 X PP G, Bi2,d = B/(10^'^G). 

§The magnetic fields (in unit of 10^^ G) derived from spectral features as cyclotron resonant scattering. No gravitation 
redshift is included here. 

tAn anomalous X-ray pulsar. 
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Iwasawa et al. (1992); 6, Marsden et al. (2001); 7, Torii et al. (1999); 8, Gotthelf et al. (2000). 
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Fig. 1. — The contour of the relative error of period for 1E1207. 4-5209 in a model with self- 
similar accretion. The star sign indicates the parametric position we choose for a reasonable 
neutron star and the fossil disk around it. The disk mass m is in Mq and the polar magnetic 
field B,o in IQi^G. 
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Fig. 2. — The period evolution in a model with self-similar accretion for the four pulsars 
labelled. The parameter sets used to calculate the curves are also hsted. The solid part 
of the curves indicates that the accretion torque works (i.e., the magnetospheric radius is 
smaller than that of light cyhnder, < tl), while the dashed part of PSR B 1757-24 means 
that the pulsar and the fossil disk could evolve independently (i.e., > tl). 



